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SUMMARY 


A simple method for extracting the period and damping and the ratios 
of all variables of the Dutch roll motion of an airplane is found by 
arranging the lateral equations of motion in such form and order that a 
rapidly convergent iterative solution may be obtained. The method is 
proposed in order to circumvent the necessity for the solution of the 
classic biquadratic characteristic equation. Because of the simplicity 
of this procedure , the iterative method is believed to be particularly 
useful when no extensive computing facilities are available, though it 
may be used to reduce computation time on any type of digital computing 
equipment. 

Primary effects of variation of the important stability derivatives 
on the period and damping can be seen more clearly from the iterative 
method than from the fourth-order characteristic equation. 


INTRODUCTION 


The importance of the short-period lateral oscillation or Dutch roll 
to the handling qualities of airplanes is recognized. (See, for example, 
refs. 1 and 2.) Although the complete solution of the lateral equations 
of motion to obtain transient response has long been possible by classi- 
cal methods, the work involved was so great that most analytical studies 
were devoted to describing the period and damping and defining the 
stability boundaries of this oscillation (refs. 3 and 4). Through the 
use of Laplace transformations, systematized, though still rather tedious, 
solutions for the complete transient motion to a given disturbance have 
been presented in references 5 and 6. 

Pilots have indicated (ref. 7) that Dutch roll characteristics are 
adequately described if the period and damping and the amplitude of the 
roll-to-sideslip ratios are known. The purpose of the present paper is to 
present a simple iterative method for just such a description of the Dutch 
roll that circumvents the necessity for the solution of the classic biquad- 
ratic characteristic equation. 
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Computational procedures and uses of the iterative method are dis- 
cussed. Numerical solutions for Dutch roll characteristics of representa- 
tive airplanes are given. 
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yawing angular velocity, ch|r/dt, radians/sec 



4 

s 

s 

t 

V 

v 

W 

X,Y,Z 

P 

T) 


P 


NACA TN 3754 

wing area, sq ft 

nondimens ional time parameter based on span, Vt/b 
time, see 
airspeed, ft/sec 

sideslip velocity about lateral axis, ft/sec 
weight of airplane, lb 

stability coordinate axes (defined in fig. 1) 

angle of sideslip, sin”- 1 - •£, radians 

inclination of principal longitudinal axis of airplane 
with respect to flight path, positive when principal 
axis is above flight path at nose, deg 

relative density factor, m/pSb 

mass density of air, slugs/cu ft 

angle of roll, radians 

angle of yaw, radians 

DEVELOPMENT OF METHOD 
Equations 


A simple method for extracting the period and damping and the ratios 
of all variables of the Dutch roll oscillation is found by arranging the 
lateral equations of motion in such form and order that a rapidly con- 
vergent iterative solution may be obtained. In this section, the original 
equations of motion are given in standard and then in modified form and 
are operated on to yield the iterative equations expressed as 

(1) The ratio of roll to yaw 0/\[r as a function of the operator D 

(2) The ratio of sideslip to yaw p/t as a function of 0/\jr and D 

(3) A quadratic in the operator D as a function of 0/\jf and p/i|r 
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Assumptions . - The usual assumptions of lateral-stability theory 
are made : 

(1) There are only three degrees of freedom: sideslip p, roll 0, 

and yaw \|r* 

(2) All aerodynamic controls are fixed. 

(3) Small classic perturbations are allowed. 

Modifications of standard equations of lateral motion .- The linear, 
second-order, simultaneous differential equations of lateral motion 
referred to stability axes (see fig, 1) for the condition of controls 
fixed are as follows: 

Side force: 

p(Wd - Cy p ) + 0 (-C L - |c Yp D) + 2pD - iCy r D) = 0 (l) 


Rolling moment: 

P ("%) + " 2% D ) + H-W 2 " 2° Z r D ) = 0 (2 > 


Yawing moment : 

P (~ Cn p) + ^ ( _2 ^ K XZ r)2 ” + ^ ( 2 ^ k Zs 2e)2 “ l^r 0 ) = ° ^ 

In the process of solving these equations by classical methods (ref. 3)> 
the characteristic equation is obtained by setting the determinant of the 
coefficients of equations (l) to (3) equal to zero. The characteristic 
equation is a biquadratic in the operator D, and it has been shown 
(ref. A) that the roots of the characteristic equation describe the 
damping of any motion present and the period if the motion is oscillatory. 

The operator D may be considered a variable in the three standard 
lateral equations of motion and still represent damping and the period if 
the motion is oscillatory. On this basis, there are now four variables: 

P, 0, i|r, and D. However, the operator D in equations (l) to (3) 
never appears except as a coefficient of one of the other variables. If 
each of the equations be divided by one of the freedoms (e.g., the 
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number of unknowns is again reduced to three: p/i|r, 0/t> and D. 

Dividing equations (l), (2), and (3) each by i|r yields the following 
equations : 

Side force: 

- C Y? ) + |(-C L - |c Yp d) + 2|i D - |c Yr D - 0 (k) 

Rolling moment: 

£('%) * - £%!>) - a^xzo 2 - = ° (5) 

Yawing moment: 

t (“S) + * (-^XZ 02 - 2 °”i. D ) + - I S: D ’ ° (6) 


Equations (4) to ( 6 ) contain enough information to define adequately the 
Dutch roll. (See ref. 7*) In this form, however, solution is admittedly 
difficult. 

Iterative equations .- The modified equations of lateral motion 
(eqs. (4) to( 6 )) will be operated on to yield the Iterative equations. 

The first iterative equation, the ratio of roll to yaw 0/ty, may be 
found by solving the roll and yaw equations (eqs. ( 5 ) and ( 6 )) simul- 
tane ous ly . Thus , 


i . S - C »B + %«*■*) (7) 

t ’ -S i% + % 1% + 2 ^(V Xs2 + % Kxz ) 

It may be noted that equation ( 7 ) is the simplest expression for any 
ratio formed from 3, 0, and ij/ inasmuch as it is the only expression 

in the lateral system that contains only first-order terms in the 
operator D and constant terms. . 
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The ratio of sideslip to yaw p/\|r may he found by rearrangement 
of terms of the remaining lateral equation, the side-force equation 
(eq. (4)). The ratio j3/\jr expressed in terms of the ratio 0/\|r, the 
operator D, and constant terms is 

3 fa - ¥**)» + l fa^ - 4%) 

-2pD + 


The remaining step is to write the quadratic in 
of 0/\|r and (3/ty, where D is the Dutch roll root, 
yawing and rolling equations contribute significantly to the Dutch roll 
oscillation, writing both the rolling -moment and the yawing-moment 
equations (eqs. ( 5 ) and (6)) as quadratics in the operator D and 
adding yields the desired quadratic in D: 


( 8 ) 


D as a function 
Since both the 


D 2 ( 2| iK Zs 2 K Xs 2 
(| c n p K Xs 2 + 





+ C 


2 " 2° Z r KxZ ” 

ifx z)f - 0 ' 


(9) 


It should be noted that before the addition the rolling -moment equation 
was multiplied by and the yawing-moment equation was multiplied by 

K Xs 2 . This procedure simplified considerably the resulting expression 
because it eliminated the complex parameter 0/t as a factor of D 2 that 
appeared in equations ( 5 ) and ( 6 ). In a somewhat arbitrary fashion, this 
procedure also actually "weighted" the importance of the contributions of 
the yawing and rolling equations to the Dutch roll roots. In the limiting 
case of zero principal-axis tilt, the weighting factor K-j^ is equal to 

zero and equation ( 9 ) reduces simply to the yawing -moment equation, equa- 
tion (6). Equation ( 9 ) also reduces to 




D = i 


( 10 ) 
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for two degrees of freedom (yaw and roll) where the system is undamped, 

3 = -4, and second-order terms are neglected. Equation (10) is identical 
to the undamped natural frequency for two degrees of freedom as defined 
in reference 8 and thus is further justification for the weighting chosen. 


Method of Solution 


The three iterative equations derived from the lateral equations of 
motion are 0/i]r in terms of D (eq. (7))> p/t in terms of 0/\|r and 
D (eq. (8)), and D in terms of 0/ijr and 3/ty (eq. (9)) wherein (/[if, 
3/i|r, and D are all complex quantities. 

Computational procedure .- The ratio /(if is estimated from equa- 
tion ( 7 ) with an assumed value of D. The ratio p/\|f is then determined 
from equation (8) with the assumed value of D and the estimated value of 
(ft /if from equation ( 7 ) being used. Next, a new value of D is determined 
from equation ( 9 ) through use of the estimated values of 0/\jr from equa- 
tion ( 7 ) and p/\|r from equation (8). This new value of D is used in 
equation ( 7 ) in place of the original assumed value, and thus the itera- 
tive process is well started. 


A suggested original assumption for the value of D is the undamped 
period for the single degree of freedom in yaw wherein the airplane center 
of gravity travels in a straight line. This well-known approximation is 



It has been found through experience that the amplitude and phase angle 
of the roll-to-yaw ratio /[if appear to be negligibly affected by 
damping and axe insensitive to the period of the Dutch roll above the 
value of D given by equation (ll). Thus, the actual choice of the 
initial D at the start of the iterative process is unimportant. Any 
reasonable approximation that presupposes an oscillation should produce a 
solution. 

The Dutch roll roots of more than 70 airplanes were calculated by 
the iterative method and checked by the conventional method. Table I 
gives some idea of the range of airplane variables investigated. As a 
result of the experience gained from these calculations, answers of 
engineering usefulness are found with the iterative method by merely 
substituting an approximate value of D in equations ( 7 ) and (8) and 
then solving equations ( 7 ) to ( 9 ) for //if, 3/t> and D, respectively. 
However, it is recommended that the process be iterated at least once 
because, if the iterated values of //if, 3/ty, and D agree with the 
values that were first calculated, the work is obviously numerically 
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correct. Two possibilities exist if-Sthe values are not in close agree- 
ment - that of numerical error and that of failure to converge. 

Infrequently, convergence does not occur (e.g., static directionally 
unstable airplanes and light airplanes that have extreme damping in yaw), 
and the conventional method must be used to evaluate the variables. 

The equations of the iterative method have been derived from the 
conventional equations of motion without any additional assumptions. As 
a result, it is impossible for the iterative method to converge on 
answers different from those calculated by the conventional method. 

Uses for the method .- The proposed iterative method is designed to 
circumvent the necessity for the solution of the biquadratic characteris- 
tic equation. The biquadratic is in effect replaced by a quadratic with 
complex coefficients (eq. ( 9 )) that may be easily iterated. Thus, the 
iterative method is somewhat simpler mathematically and requires only 
about one -third the computational time of the conventional method. In 
addition, evaluation of the ratios of roll to yaw and sideslip to yaw 
which are necessary to describe adequately the Dutch roll is inherent 
with the iterative method. 

Equation ( 9 ) may be rewritten as 

D 2 + D " I^Xs 2 " - (^n/xs 2 + | C Z p K Xz)f' 

2p(k Z s 2k Xs 2 - K^ 2 ) 


( C n p %3 2 + C Z p K Xz)f 

^As 2 “ ^Z 2 ) 


( 12 ) 


Equation (12) is an extremely interesting equation in that it contains all 
the terms in a three-degree -of-freedom system that affect the Dutch roll, 
expressed as a quadratic in the operator D. Primary effects of variation 
of the important stability derivatives on the period and damping can be 
seen more clearly from equation (12) than from the fourth-order character- 
istic equation. 
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REPRESENTATIVE SOLUTIONS 


Table II presents the characteristics of four representative air- 
planes. These airplanes were examined by means of the simple iterative 
method to determine their Dutch roll characteristics. 

Example solutions by the iterative method are given in table III. 

A conventional modern bomber, an extreme -altitude fighter, and a sonic 
interceptor represent typical solutions for present and proposed air- 
planes, and a hypothetical delta-wing light airplane in the landing con- 
dition represents an atypical solution in which the iterative method 
does not converge. The high rate of convergence of the typical 
solutions of the Dutch roll characteristics is apparent. The atypical 
case illustrates the ample warning (lack of convergence) that is present 
when the iterative solution should not be used. 


CONCLUSIONS 


A simple method for calculating the characteristics of the Dutch 
roll motion of an airplane has been obtained by arranging the lateral 
equations of motion in such form and order that an iterative process is 
quickly convergent. The iterative method is believed to be particularly 
useful when no extensive computing facilities are available, although it 
may be used to reduce computational time on any type of digital computing 
equipment. Experience gained from the calculation of the Dutch roll 
characteristics of more than JO airplanes by the Iterative method as 
compared with the conventional method has indicated that 

1. About one-third the computational time Is required for the 
iterative method as was required for the conventional method. 

2. The arithmetical processes are simpler in the iterative method 
because the need for the solution of a biquadratic equation is avoided. 

3 . In the iterative method the arithmetical processes are iterative, 
which permit running numerical checks. 

4. Evaluation of the ratios of roll to yaw and sideslip to yaw is 
inherent in the iterative process. 
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5- Primary effects of aerodynamic derivatives on the Dutch roll 
roots can he more readily seen by use of the iterative method. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., June 8, 195 6 * 
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TABLE I 


RANGE OF AIRPLANE VARIABLES CONSIDERED FOR WHICH 
SATISFACTORY SOLUTIONS WERE OBTAINED 
BY MEANS OF ITERATIVE METHOD 


Primary variable 

Numerical range 

^s 2 

0.2 to 1.0 



P 

12 to 275 

T], deg 

-5 to 5 

% 

-1.2 to 3.6 

% 

-O.7O to 0.90 

c z 

b T 

-2.90 to 3.10 


C z 

-7.3O to 0.20 


% 

0.03 to 0.25 

c z 

Z j 3 

-0.22 to 0 




a t= 
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L * * * ' 

h , ft . . 
V, ft/sec 
ft . . 


Airplane 

Bomber 

Extreme - 
altitude 
fighter 

Sonic 

interceptor 

Hypothetical 
delta-wing 
airplane in 
landing 
condition 

0.0311 

O.OI56 

0.010 

0.030 

0.072 

O.I56 

O.065 

O.0679 

0 

0.0020 

-0.004 

-0.030 

31.83 

182 

65A 

II.85 

0.443 

0.49 

0.323 

1.0 

35,000 

50.000 

32,000 

0 

700 

776 

691 

99-8 

116 

25 

36.6 

38.1 

-0.61 

-O.58 

-0.28 

-0.286 

0 

0 

0 

0 

0 

0 

0 

0 

-0.14 

-0.18 

-0.08 

-0.0573 

-0.44 

-0.33 

-0.15 

-0.02 

0.149 

0.23 

0.40 

0.60 

0.12 

0.25 

0.028 

0.0573 

-0.0276 

-0.050 

-0.20 

-0.20 

-0.156 

-0.069 

-0.40 

-1.10 








NACA TN 3754 


TABLE III 


REPRESENTATIVE SOLUTIONS 


Assumptions 





Airplane 


Conventional bomber 


4 >h 

p/+ 

D 

Iterated values 

-1.755 - 1-2551 
-1.840 - 1.2351 
-1.825 - 1.2431 

-1.053 + 0.013171 
-1.053 + 0.01771 
-1.053 + 0 . 016921 

-0.00491 + 0.16791 
-0.00441 + 0.16791 
-0.00452 + 0.16791 

Exact values 

-1.825 - 1.2421 

-1.053 + 0. 016921 

-0.00447 + 0.16791 


Extreme -altitude fighter 


ih 

p h 

D 

Iterated values 

-6.23 - 2.951 
-6.17 - 2.841 
-6.17 - 2.851 

-1.062 + 0.10081 
-1.064 + 0.09681 
-1.064 + 0.09701 

0.00271 + 0.06651 
0.00256 + 0.06651 
0.00257 + 0.06651 

Exact values 

-6.17 - 2.851 

-1.064 + 0.09701 

0.00258 + O.O665I 


Sonic interceptor 



P h 

D 

Iterated values 

-2.58 - 2.271 
-4.4l - 2.541 
-4.33 - 2.461 

-1.102 + 0.07001 
-1.045 + 0.07751 
-1.047 + 0.07491 

-0.01395 + 0.11981 
-O.OO852 + 0.11821 
-0.00886 + 0.11811 

Exact values 

-4.33 - 2.46i 

-1.047 + 0.07511 

-O.OO885 + 0.11811 


Hypothetical delta-wing light airplane in landing condition 


0/+ 

p/t 

D 

Iterated values 


-1.724 + 0.2251 
-1.411 + 0.2311 
-1.039 + 0.2341 

-O.97O + 0.323i 
-O.781 + 0.08621 
-0.764 - 0.02271 

-0.1845 

-0.1498 

-0.1389 

+ 0.11231 
+ 0.02861 
+ 0.0006171 

Exact values 


1.722 + 0.589i I -0.451 + 0.3851 | -0.0647 + 0.0822i 

























































